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Using data samples collected at center-of-mass energies of y/s = 4.009, 4.226, 4.257, 4.358, 4.416, 
and 4.599 GeV with the BESIII detector operating at the BEPCII storage ring, we search for the 
isospin violating decay E(4260) —> J/ifrin 0 . No signal is observed, and upper limits on the cross 
section a(e + e~ —> J/ipr/n 0 ) at the 90% confidence level are determined to be 3.6, 1.7, 2.4, 1.4, 0.9, 
and 1.9 pb, respectively. 

PACS numbers: 14.40.Rt, 13.66.Be, 14.40.Pq, 13.20.Gd 


I. INTRODUCTION 


to open charm decay modes |2j. y(4260) is a vector 


The y(4260) charmoniumlike state was first observed 
in its decay to n + TT~J/’tl> (Tj and has a small coupling 
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( J PC = 1 ) state that is only barely observable as an s- 

channel resonance in e + e _ collisions and that appears at 
an energy where no conventional charmonium state is ex¬ 
pected. Since its discovery, many theoretical studies have 
been carried out considering the Y (4260) as a tetraquark 
state Q, D\D or DqD* hadronic molecule [IJ, hybrid 
charmonium Q, baryonium state etc. 

Recently, in the study of Y(4260) —> J/ip, a 

charged charmoniumlike structure, the Z c ( SOOO)^, was 
observed in the i invariant mass spectrum by 
the BESIII 3 and Belle experimentsand confirmed 
shortly thereafter with CLEO-c data [9]. In the molecule 
model M , the F(4260) is proposed to have a large D\D 
component, while Z c (39 00) ± has a DD* component. 

BESIII recently reported the observation of e + e _ —>■ 
7^(3872) —y 77 r + 7 r~ J/i/ [TT|. The cross section mea¬ 
surements strongly support the existence of the radiative 
transition V(4260) —> 7X(3872). One significant feature 
of the X(3872) that differs from conventional charmo¬ 
nium is that the decay branching fraction of _X_[3872) to 
7 r + 7 r - 7 T° J/'0 is comparable to n + n~J/t/j [l2j, Ol, so the 
isospin violating process occurs on a large scale. 

Isospin violating decays can be used to probe the 
nature of heavy quarkonium. The hadro-charmonium 
model [l4|] and tetraquark models Ei predict that 
the reaction T(55) —► rjTr°+ bottomonium should be ob¬ 
servable. The tetraquark model E2 also predicts that Z ° 
can be produced in Y (4260) —> J/^-k 0 with Z° decaying 
into J/V>7T° and possibly J/ipV in the presence of sizable 
isospin violation. The molecular model |l 8 j predicts a 
peak in the cross section of V(4260) -7 J/ipriit 0 at the 
D\D threshold and a narrow peak in the J/^uy invariant 
mass spectrum at the DD* threshold. 

In this paper, we present results on a search for the 
isospin violating decay V(4260) —> J/i/v7 r°, with J/tj) -7 
e + e~ /p + p,~, 7T° -7 77 , and p —> 77 (the other decay 
modes of ry are not used due to much lower detection effi¬ 
ciency and branching fraction), based on e + e _ annihila¬ 
tion data collected with the BESIII detector operating at 
the BEPCII storage ring [l9| at center-of-mass energies 
of Vs = 4.009, 4.226, 4.257, 4.358, 4.416, and 4.599 GeV. 


II. BESIII DETECTOR AND MONTE CARLO 
SIMULATION 


The BESIII detector, described in detail in Ref. TT)!. 
has a geometrical acceptance of 93% of 47 r. A small¬ 
cell helium-based main drift chamber (MDC) provides 
a charged particle momentum resolution of 0.5% at 1 
GeV/c in a 1 T magnetic field and supplies energy- 
loss ( dE/dx ) measurements with a resolution of 6 % for 
minimum-ionizing pions. The electromagnetic calorime¬ 
ter (EMC) measures photon energies with a resolution of 
2.5% (5%) at 1.0 GeV in the barrel (end caps). Particle 
identification is provided by a time-of-flight system with 


a time resolution of 80 ps (110 ps) for the barrel (end 
caps). The muon system (MUC), located in the iron flux 
return yoke of the magnet, provides 2 cm position reso¬ 
lution and detects muon tracks with momentum greater 
than 0.5 GeV/c. 


The GEANT4-based [2Cj Monte Carlo (MC) simulation 
software BOOST [21] includes the geometric description 
of the BESIII detector and a simulation of the detector 
response. It is used to optimize event selection criteria, 
estimate backgrounds, and evaluate the detection effi¬ 
ciency. For each energy point, we generate large signal 
MC samples of e + e _ -7 J/i/pn 0 , J/ip -7 e + e~/p, + p~, 
V —■ > 77 , and n° —> 77 uniformly in phase space. Ef¬ 
fects of initial state radiation (ISR) are simulated with 


KKMC 


where the Born cross section of e + e 


J/i/p 7 T° is assumed to follow a V(4260) Breit—Wigner 
line shape with resonance parameters taken from the Par¬ 
ticle Data Group (PDG) [23|. Final state radiation ef¬ 
fects associated with charged particles are handled with 
PHOTOS [2I. 


To study possible backgrounds, a MC sample of inclu¬ 
sive Y (4260) decays, equivalent to an integrated luminos¬ 
ity of 825.6 pb -1 , is also generated at V s = 4.260 GeV. In 
these simulations, the Y (4260) is allowed to decay gener- 
ically, with the main known decay channels being gener¬ 
ated using EVTGEN 12511 with branching fractions set to 
world average values |23|. The remaining events associ¬ 
ated with charmonium decays are generated with LUND- 
CHARM [26|, while continuum hadronic events are gen¬ 
erated with PYTHIA [27j |. QED events (e + e _ -7 e + eN 
/i + yr - , and 77 ) are generated with KKMC [22]. Back¬ 
grounds at other energy points are expected to be simi¬ 
lar. 


III. EVENT SELECTION 


Events with two charged tracks with a net charge of 
zero are selected. For each good charged track, the polar 
angle in the MDC must satisfy |cos0| < 0.93, and the 
point of closest approach to the e + e _ interaction point 
must be within ±10 cm in the beam direction and within 
±1 cm in the plane perpendicular to the beam direction. 
The momenta of leptons from the J/i[) decays in the lab¬ 
oratory frame are required to be larger than 1.0 GeV/c. 
E/p is used to separate electrons from muons, where E 
is the energy deposited in the EMC and p is the mo¬ 
mentum measured by the MDC. For electron candidates, 
E/p should be larger than 0.7, while for muons, it should 
be less than 0.3. To suppress background from events 
with pion tracks in the final state, at least one of the two 
muons is required to have at least five layers with valid 
hits in the MUC. 

Showers identified as photon candidates must satisfy 
fiducial and shower quality as well as timing require¬ 
ments. The minimum EMC energy is 25 MeV for barrel 
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showers (| cos0| < 0.80) and 50 MeV for end cap showers 
(0.86 < |cos0| < 0.92). To eliminate showers produced 
by charged particles, a photon must be separated by at 
least 5 deg from any charged track. The time information 
from the EMC is also used to suppress electronic noise 
and energy deposits unrelated to the event. At least four 
good photon candidates in each event are required. 

To improve the momentum resolution and reduce the 
background, the event is subjected to a four-constraint 
(4C) kinematic fit under the hypothesis e + e _ —> 
7777 l + l~ (l = e/p,), and the y 2 is required to be less 
than 40. For events with more than four photons, the 
four photons with the smallest y 2 from the 4C fit are 
assigned as the photons from 77 and n°. 

After selecting the 7777 l + l~ candidate, scatter plots 
of M(77) with all six combinations of photon pairs for 
events in the J/ip signal region (3.067 < M(l + l~ ) <3.127 
GeV/c 2 ) for data at yfs = 4.226 and 4.257 GeV are shown 
in the left two panels of Fig. |T] Distributions of M ( l + l ~) 
for events in the tt 0 tt 0 signal region (both photon pairs 
satisfy |M(y7) — m n o| < 10 MeV/c 2 ) for data at yfs = 
4.226 and 4.257 GeV are shown in the right two panels of 
Fig. [U Clear J/'p peaks are observed, corresponding to 
7r°7r° J/ip events. To remove this 7r°7r° J/ip background, 
events with any combination of photon pairs in the 7 r° 7 r° 
region of the scatter plot are rejected. 



M(yy) (GeV/c 2 ) M(IT) (GeV/c 2 ) 


FIG. 1 : Scatter plot of Mfy 7) with all six combinations 
for events in the J/ip signal region (left) and distribution 
of for events in the 7 r° 7 r° signal region (right) for 

data at y/s = 4.226 GeV (top) and 4.257 GeV (bottom). 

After rejecting the 7r°7r°J/^> background, we choose 
the combination of photon pairs closest to the rjn° signal 

region by minimizing ^ p'^’G' 1 ^ 3 ) —^-2.| 2 -|-|- 1 13 -| 2 , 

where <j v and a^o are the 77 and 7r° resolutions obtained 
from the signal MC, respectively. The scatter plots of 
M (77) with the combination closest to the r/n 0 signal re¬ 
gion for events in the J/ip signal region for data at yfs = 
4.226 and 4.257 GeV are shown in the top two panels of 


Fig. [ 2 ] No cluster of r/n 0 events is observed in the J/^ sig¬ 
nal region, with a vertical band for n° —> 77 clearly visi¬ 
ble, but no prominent band for rj —> 77 is observed. The 
projections of the scatter plots on M (7172) with M (7374) 
in the 7r° signal region (|-M/7374) — m, r o | <10 MeV/c 2 ) 
and projections on M(7374) with -M/7172) in the 77 signal 
region [\M{^1^-2) — mf < 30 MeV/c 2 ) for data are shown 
in the middle and bottom panels of Fig. ( 2 j respectively. 



FIG. 2: Scatter plot of .M/ 77 ) for the combination closest 
to the rjir 0 signal region for events in the J/t/> signal region 
(top), projection of the scatter plot on M( 7 i 7 2 ) with 
1 ^( 7374 ) in 7T° signal region (middle), and projection of 
the scatter plot on 114 ( 7374 ) with M( 7172 ) in 77 signal 
region (bottom) for data at yfs = 4.226 GeV (left) and 
4.257 GeV (right). 

The background for e + e _ — > J/firiTT 0 is studied us¬ 
ing the inclusive MC sample at yfs = 4.260 GeV. Af¬ 
ter imposing all event selection requirements, there are 
two background events from e + e _ — > 7r°7r° J/if) and nine 
background events arising from e + e - — > 7isr'0 , ) 7 isrV’ ,/ ; 
and 7 isrt/>( 4040). No other background survives. The 
background can be evaluated with r/ir° sideband events. 
Distributions of M(l + l~) for events in the tjtt 0 signal re¬ 
gion for data at yfs = 4.226 and 4.257 GeV are shown in 
Fig. [3] Distributions of M(l + l~) for events correspond¬ 
ing to the normalized two-dimensional tjtt 0 sidebands 
are shown as shaded histograms. The 77 sideband re¬ 
gions are defined as 0.3978 < M{ 7172 ) < 0.4578 GeV/c 2 
and 0.6378 < M( 7 172 ) < 0.6978 GeV/c 2 . The 7 r° side¬ 
band regions are defined as 0.0849 < M( 7374 ) < 0.1049 
GeV/c 2 and 0.1649 < M( 7 374 ) < 0.1849 GeV/c 2 . The 
counted number of observed events in the J/iI>t]tt 0 signal 
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region N° bs and number of background events estimated 
from rjTT° sidebands 7V bkg are listed in Table E 
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FIG. 3: Distributions of M(l + l ~) for events in r/n° sig¬ 
nal region and sideband regions for data at y/s = 4.226 
GeV (left) and 4.257 GeV (right). The error bars are the 
M(l + l ~) distributions for events in the 7771 -° signal region, 
and the shaded histograms are those in the t]tt° sideband 
regions. 


IV. CROSS SECTION UPPER LIMITS 


Since no signal above the background is ob¬ 

served, upper limits on the Born cross section of e + e _ —> 
J/^r? 7 r° at the 90% C.L. are determined using the for¬ 
mula 

N up 

" C(l + 6 r )(l + 5 v ){e ee B ee + ewBw)B*°&> ’ 

where N/ P seived is the upper limit on the number of sig¬ 
nal events; C is the integrated luminosity; (1 + 5 r ) is the 
radiative correction factor, which is taken from a QED 
calculation assuming the e + e~ —> J/’/rj'K 0 cross section 
is described by a Y (4260) Breit—Wigner line shape with 
parameters taken from the PDG [HI; (1 Td 11 ) is the vac¬ 
uum polarization factor including leptonic and hadronic 
parts and taken from a QED calculation with an accuracy 
of 0.5% j28|; e ee and are the efficiencies for e + e~ and 
/i + ^ _ modes, respectively; B ee and B w are the branch¬ 
ing fractions of J/if) —y e + e _ and J// —> [HI, 

respectively; and B n and B v are the branching fractions 
of r] —> 77 and ir° —> 77 [HI, respectively. 

The efficiency corrected upper limit on the number of 

jV up i 

signal events N up = is estimated with N° as 

and iV bkg using the profile likelihood method, which is 
implemented by TRolke in the ROOT framework [2f| • The 
calculation for obtaining N up includes the background 
fluctuation and the systematic uncertainty of the cross 
section measurement. The background fluctuation is as¬ 
sumed to follow a Poisson distribution. The systematic 
uncertainty of the cross section is taken as a Gaussian 
uncertainty. 

The systematic uncertainty of the cross section mea¬ 
surement in Eq. m includes the luminosity mea¬ 
surement, detection efficiency, and intermediate decay 
branching fractions. The systematic uncertainties of the 


luminosity, track reconstruction, and photon detection 
are 1.0% yd]], 1.0% per track (30|, and 1.0% per pho¬ 
ton [3 1 |, respectively. The systematic uncertainties from 
the branching fraction of 7r° and 77 decays are taken 
from the PDG (23j. These sources of systematic uncer¬ 
tainty, which are summarized in the top part of Table QT] 
are common for e + e _ and modes. The follow¬ 

ing sources of systematic uncertainty, which are uncor¬ 
related for the e + e _ and fi + fi~ modes, are summarized 
in the bottom part of Table [0] The systematic uncer¬ 
tainty from the branching fraction of J/ip decay is taken 
from the PDG [HI- The systematic uncertainty from 
the requirement on the number of MUC hits is 3.6% and 
estimated by comparing the efficiency of the MUC re¬ 
quirement between data and MC in the control sample 
e + e _ —> n 0 TT 0 J/ilj at y/s = 4.257 GeV. The systematic 
uncertainty from the requirement of the J/ip signal region 
is estimated by smearing the invariant mass of l + l~ of the 
signal MC with a Gaussian function to compensate for 
the resolution difference between the data and MC when 
calculating the efficiency. The parameters for smearing 
are determined by fitting the J/ip distribution of data 
with the MC shape convoluted with a Gaussian function 
for the control sample e + e _ —> n°n°J/i/j. The difference 
in the detection efficiency between signal MC samples 
with and without the smearing is taken as the system¬ 
atic uncertainty. The systematic uncertainty from the 
MC model is estimated by generating a MC sample with 
the angular distribution of leptons determined from the 
J/ip data. The systematic uncertainty due to kine¬ 
matic fitting is estimated by correcting the helix parame¬ 
ters of charged tracks according the method described in 
Ref. [32|, where the correction factors are obtained from 
the control sample if)' —> JXcJ and the difference in the 
detection efficiency between with and without making 
the correction to the MC is taken as the systematic un¬ 
certainty. The uncorrelated systematic uncertainties for 
the electron and muon channels are combined by taking 
the weighted average with weights e ee B ee and e^B^, re¬ 
spectively. The total systematic uncertainty is obtained 
by summing all the sources of the systematic uncertainty 
in quadrature. 

The systematic uncertainty on the size of the back¬ 
ground is estimated by evaluating N up with different sig¬ 
nal and sideband regions for 77 and 7 r°. The most con¬ 
servative N up is taken as the final result, as listed in 
Table E The upper limits on the Born cross section of 
e + e _ —> J/ipijK 0 (ctJjl™) assuming it follows a V(4260) 
Breit—Wigner line shape are listed in Table E 

For comparison, the radiative correction factor and de¬ 
tection efficiency have been recalculated assuming the 
e + e _ —> J/^rjn 0 cross section follows alternative line 
shapes. If the cross section follows the line shape of the 
V(4040), the upper limit on the Born cross section is 4.1 
pb at y/s = 4.009 GeV. For a V(4360) line shape, it is 
1.6 pb at y/s = 4.358 GeV. For a V(4415) line shape, it 
is 1.5 pb at y/s = 4.358 GeV and 1.0 pb at y/s = 4.416 
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TABLE I: Results on e + e“ —> J/^-k 0 . Listed in the table are the integrated luminosity £, radiative correction 
factor (l+V) taken from QED calculation assuming the Y(4260) cross section follows a Breit—Wigner line shape, 
vacuum polarization factor (1+5’'), average efficiency (e ee B ee + number of observed events N ohs , number 

of estimated background events V bkg , the efficiency corrected upper limits on the number of signal events 7V up , and 
upper limits on the Born cross section er®” 11 (at the 90 % C.L.) at each energy point. 


y/i (GeV) C (pb’ 1 ) 

( 1 +V) 

(1+V) (e 

eegee _|_ y obs 

jV bk s 

N up 

(Pb) 

4.009 

482.0 

0.838 

1.044 

2.1 ± 0 . 1 (sys.) 

5 

1 

598.1 

3.6 

4.226 

1047.3 

0.844 

1.056 

2.2 ± 0 . 1 (sys.) 

12 

11 

592.9 

1.7 

4.257 

825.6 

0.847 

1.054 

2.2 ± 0 . 1 (sys.) 

12 

8 

654.1 

2.4 

4.358 

539.8 

0.942 

1.051 

2.2 ± 0 . 1 (sys.) 

5 

4 

283.2 

1.4 

4.416 

1028.9 

0.951 

1.053 

2.3 ± 0.1(sys.) 

5 

6 

342.7 

0.9 

4.599 

566.9 

0.965 

1.055 

2.4 ± 0.1(sys.) 

6 

3 

418.4 

1.9 


TABLE II: Systematic uncertainties in the J/ipT]TT° cross section measurement at each energy point (in %). The 
items in parentheses in the bottom part of the table are the uncorrelated systematic uncertainties for the e + e - 
(first) and (second) modes. 


Sources/ y/s (GeV) 

4.009 

4.226 

4.257 

4.358 

4.416 

4.599 

Luminosity 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

MDC tracking 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

Photon reconstruction 

4.0 

4.0 

4.0 

4.0 

4.0 

4.0 

B( 7T° -A 77 ), B (7 — > 77 ) 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

B{j/ip -a i+r) 

(0.5, 0.5) 

(0.5, 0.5) 

(0.5, 0.5) 

(0.5, 0.5) 

(0.5, 0.5) 

(0.5, 0.5) 

MUC hits 

(0, 3.6) 

(0, 3.6) 

(0, 3.6) 

(0, 3.6) 

(0, 3.6) 

(0, 3.6) 

J/ip mass resolution 

(0.2, 1.3) 

( 0 . 8 , 1 . 2 ) 

(0.5, 1.3) 

(0.2, 0.7) 

(0.7, 1.6) 

( 0 . 1 , 0 . 6 ) 

Decay model 

(1.5, 1.9) 

(0.9, 1.1) 

(0.4, 0.6) 

(0.2, 0.7) 

(0.7, 0.2) 

( 0 . 2 , 0 . 2 ) 

Kinematic fitting 

(1.2, 0.9) 

( 1 . 1 , 1 . 2 ) 

(0.9, 0.9) 

(0.7, 1.2) 

( 1 . 1 , 1 . 0 ) 

(1.0, 1.4) 

Total 

5.3 

5.3 

5.2 

5.2 

5.3 

5.2 


GeV. For a Y (4660) line shape, it is 2.0 pb at y/s = 4.599 
GeV. 

It is also possible to set upper limits on e + e“ —> 
Z°7r° -A J/ifjri'K 0 . The number of observed events and 
number of estimated background events in the Z® signal 
region (3.850 < M(J/ipr]) < 3.940 GeV/c 2 ) are 7 and 4 
±2, respectively, at y/s = 4.226 GeV, and 8 and 3 ± 2, 
respectively, at y/s = 4.257 GeV. The upper limit on 
cr(e + e~ -A Z®n° -A J/i/jr/n 0 ) is determined to be 1.3 pb 
at y/s = 4.226 GeV and 2.0 pb at y/s = 4.257 GeV, where 
only the statistical uncertainty is given. Compared to the 
measured cross section of e + e _ —> Z°7r° -a J/ipTr°TT° pSSij ]. 
the upper limit on the ratio of the branching fraction 

90% confidence level is 0.15 at y/s = 
4.226 GeV and 0.65 at y/s = 4.257 GeV. 

V. SUMMARY 

In summary, using data collected with the BESIII de¬ 
tector, a search for the isospin violating decay Y (4260) —1 
J/ipr]i: Q is performed. No statistically significant signal 
is observed. The Born cross sections of e + e _ -A 
at the 90% confidence level limits at y/s = 4.009, 4.226, 
4.257, 4.358, 4.416, and 4.599 GeV are determined to be 
3.6, 1.7, 2.4, 1.4, 0.9, and 1.9 pb, respectively. The up¬ 
per limits are well above the prediction for the molecule 


model [3. 
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